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ABSTRACT 


Local  and  average  heat  transfer  coefficients  were  deter- 
mined for  a  right  circular  cylinder  in  an  oscillating  flow. 
Spanwise  platinum  heater  strips  were  used  to  heat  the  cylinder 
isothermally  over  the  lower  180  degrees  from  the  front  to  rear 
stagnation  point.   The  four  inch  diameter  cylinder  was  posi- 
tioned both  normal  to  and  at  45  degrees  to  the  flow  direction. 

Data  was  gathered  for  diameter  Reynolds  numbers  from 
100,000  to  300,000.   Large  amplitude  oscillations  were  imposed 
upon  the  mean  flow  using  a  rotating  shutter  arrangement.   Fre- 
quencies of  0,  5,  10,  22,  50,  100  and  126  Hertz  were  investigated 

For  normal  flows,  local  heat  transfer  coefficients  in  the 
wake  and  average  Nusselt  numbers  were  enhanced  above  values 
for  steady  flow,  for  runs  having  large  amplitude.   In  the  45° 
flows,  no  significant  change  in  heat  transfer  was  noted  with 
oscillating  flow. 


TABLE  OF  CONTENTS 

I.  INTRODUCTION 12 

A.   NORMAL  CYLINDER 12 

II.  DESCRIPTION  OF  APPARATUS 15 

A.  MODEL 15 

B.  WIND  TUNNEL 16 

III.  EXPERIMENTAL  PROCEDURE  21 

A.  CALIBRATION  OF  MODEL 21 

B.  ADJUSTMENT  OF  SURFACE  HEATERS   22 

C.  FREESTREAM  CONDITIONS  23 

D.  DATA  REDUCTION 23 

IV.  RESULTS  AND  DISCUSSION 25 

A.  STEADY  FLOW 25 

B.  NONSTEADY  FLOW 26 

C.  NUSSELT  NUMBER  CORRELATION  28 

V.  CONCLUSIONS 29 

APPENDIX  A  -  Data  Reduction 30 

APPENDIX  B  -  Sample  Calculation  32 

APPENDIX  C  -  Radiation 36 

APPENDIX  D  -  Conduction  Losses  38 

APPENDIX  E  -  Uncertainty  Analysis  42 

APPENDIX  F  -  Predictor  -  Corrector 

Equation 48 


LIST  OF  REFERENCES  .  .  . 
INITIAL  DISTRIBUTION  LIST 


143 
146 


LIST  OF  TABLES 

I.  Summary  of  Average  Heat  Transfer  Results  103 

II.  Summary  of  Local  Heat  Transfer  Results  105 


LIST  OF  FIGURES 


Figure 
1. 
2. 
3. 

4. 
5. 

6. 
7. 
8. 

9. 
10. 
11. 
12. 
13. 

14. 

15. 


16 


17. 


18. 


Sketch  of  the  Model 50 

Cross  Sectional  Sketch  of  the  Model 51 

Photograph  of  the  Model 52 

Photograph  of  the  Wind  Tunnel  Test  Facility 53 

Photograph  of  the  Wind  Tunnel  Test  Section 

and  Instrumentation 54 

Plan  View  of  the  Wind  Tunnel 55 

Photograph  of  the  Rotating  Shutter  Valves  56 

Photograph  of  the  Model  Mounted  in  the 

Tunnel 57 

Sketch  of  Calibration  and  Working  Circuits  58 

Typical  Velocity  Waveforms  59 

Sketch  of  Teledeltos  Paper  Model  ...  60 

Sketch  of  the  Lumped  Resistance  Model  61 

Local  Heat  Transfer  Coefficients  in  Steady 

Flow 6  2 

Local  Heat  Transfer  Coefficients  in  Steady 

Flow 63 

Local  Heat  Transfer  Coefficients  in  Steady 

Flow 64 

Local  Heat  Transfer  Coefficients  in  Steady 

Flow 65 

Local  Heat  Transfer  Coefficients  in  Steady 

Flow 66 

Local  Heat  Transfer  Coefficients  in  Steady 

Flow 67 

8 


19.     Comparison  of  Published  Steady  Flow  Results 
to  Experimental  


68 


20.  Local  Heat  Transfer  Coefficients  in 
Oscillating  Flow  that  Showed  Unaltered  Results 

21.  Local  Heat  Transfer  Coefficients  in 
Oscillating  Flow  that  Showed  Unaltered  Results 

22.  Local  Heat  Transfer  Coefficients  in 
Oscillating  Flow  that  Showed  Unaltered  Results 

23.  Local  Heat  Transfer  Coefficients  in 
Oscillating  Flow  that  Showed  Unaltered  Results 

24.  Local  Heat  Transfer  Coefficients  in 
Oscillating  Flow  that  Showed  Unaltered  Results 

25.  Local  Heat  Transfer  Coefficients  in 
Oscillating  Flow  that  Showed  Unaltered  Results 

26.  Local  Heat  Transfer  Coefficients  in 
Oscillating  Flow  that  Showed  Unaltered  Results 

27.  Local  Heat  Transfer  Coefficients  in 
Oscillating  Flow  that  Showed  Unaltered  Results 

28.  Local  Heat  Transfer  Coefficients  in 
Oscillating  Flow  that  Showed  Unaltered  Results 

29.  Local  Heat  Transfer  Coefficients  in 
Oscillating  Flow  that  Showed  Unaltered  Results 

30.  Local  Heat  Transfer  Coefficients  in 
Oscillating  Flow  that  Showed  Unaltered  Results 

31.  Local  Heat  Transfer  Coefficients  in 
Oscillating  Flow  that  Showed  Unaltered  Results 

32.  Local  Heat  Transfer  Coefficients  in 
Oscillating  Flow  that  Showed  Unaltered  Results 

33.  Local  Heat  Transfer  Coefficients  in 
Oscillating  Flow  that  Showed  Unaltered  Results 

34.  Local  Heat  Transfer  Coefficients  in 
Oscillating  Flow  that  Showed  Unaltered  Results 

35.  Local  Heat  Transfer  Coefficients  in 
Oscillating  Flow  that  Showed  Unaltered  Results 


69 


70 


71 


72 


73 


74 


75 


76 


77 


79 


80 


81 


82 


.  83 


84 


36.  Local  Heat  Transfer  Coefficients  in 

Oscillating  Flow  that  Showed  Unaltered  Results   ...  85 

37.  Local  Heat  Transfer  Coefficients  in 

Oscillating  Flow  that  Showed  Unaltered  Results   ...  86 

38.  Local  Heat  Transfer  Coefficients  in 

Oscillating  Flow  that  Showed  Unaltered  Results   ...  87 

39.  Local  Heat  Transfer  Coefficients  in 

Oscillating  Flow  that  Showed  Unaltered  Results   ...  88 

40.  Local  Heat  Transfer  Coefficients  in 

Oscillating  Flow  that  Showed  Enhanced  Heat  Transfer  .  89 

41.  Local  Heat  Transfer  Coefficients  in 

Oscillating  Flow  that  Showed  Enhanced  Heat  Transfer  .  90 

42.  Local  Heat  Transfer  Coefficients  in 

Oscillating  Flow  that  Showed  Enhanced  Heat  Transfer  .  91 

43.  Local  Heat  Transfer  Coefficients  in 

Oscillating  Flow  that  Showed  Enhanced  Heat  Transfer  .  92 

44.  Local  Heat  Transfer  Coefficients  in 

Oscillating  Flow  that  Showed  Enhanced  Heat  Transfer  .  93 

45.  Local  Heat  Transfer  Coefficients  in 

Oscillating  Flow  that  Showed  Enhanced  Heat  Transfer  .  94 

46.  Local  Heat  Transfer  Coefficients  in 

Oscillating  Flow  that  Showed  Enhanced  Heat  Transfer  .  95 

47.  Local  Heat  Transfer  Coefficients  in 

Oscillating  Flow  that  Showed  Enhanced  Heat  Transfer  .  96 

48.  Local  Heat  Transfer  Coefficients  in  Oscillating 

Flow  that  Showed  Degraded  Heat  Transfer 97 

49.  Local  Heat  Transfer  Coefficients  in  Oscillating 

Flow  that  Showed  Degraded  Heat  Transfer 98 

50.  Local  Heat  Transfer  Coefficients  in  Oscillating 

Flow  that  Showed  Degraded  Heat  Transfer 99 

51.  Local  Heat  Transfer  Coefficients  in  Oscillating 

Flow  that  Showed  Degraded  Heat  Transfer 100 

52.  Variation  of  Average  Nusselt  number  with 

Reynolds  number   101 


10 


53.     Variation  of  Average  Nusselt  number  with 
Reynolds  number  


102 


11 


I .   INTRODUCTION 

A.   NORMAL  CYLINDER 

Forced  convective  heat  transfer  from  a  cylinder  in  uniform 
flow  has  been  a  topic  of  interest  to  researchers  since  the  turn 
of  the  century.   In  1914  Professor  L.V.  King  presented  his 
relationship  to  describe  the  heat  transfer  from  a  cylinder  in 
laminar  crossflow  as  a  function  of  air  velocity: 


Nu  =  A.  +  B.NB  °*5 
1     1  Re 


Although  slightly  modified  by  Collis  and  Williams  in  1959, 
it  still  serves  today  and  is  the  basic  equation  of  hot  wire 
anemometery.   Numerous  other  correlations  are  available  and 
are  a  standard  part  of  the  heat  transfer  literature  [Ref .  1] . 

A  topic  of  increasingly  greater  interest  to  researchers 
and  designers  is  the  effect  of  an  oscillating  flow  on  forced 
convection.   Gas  turbine  designers  are  intensely  interested 
in  upgrading  the  performance  of  gas  turbine  blades  through 
materials  research  and  blade  cooling.   Gas  turbine  blades 
experience  a  hot,  oscillating  flow  environment  as  they  pass 
through  the  wakes  of  upstream  blades. 

Heat  exchange  designers  involved  in  waste  heat  recovery 

from  diesel  engine  exhausts  are  also  concerned  with  the  effect 

of  flow  oscillations  on  the  exchange  of  heat  to  and  from 

cylinders . 
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Little  data  has  been  made  available  to  date  for  a  cylinder 
in  oscillating  flow.   Data  is  available  for  heat  transfer  from 
a  flat  plate  in  oscillating  flow  [Refs.  1,  2  and  3]   and  for 
an  airfoil  in  oscillating  flow  [Ref.  4].   Data  for  a  cylinder 
in  other  than  uniform  flow  is  generally  for  a  vibrating  cylinder 
with  or  without  a  net  flow  [Ref.  5]   or  for  a  cylinder  with 
screen  induced  turbulence  of  less  than  about  ten  percent  [Refs. 
6,  7,  8  and  9].   This  study  concerns  itself  with  the  effect  of 
oscillations  on  the  order  of  fifty  percent  of  the   freestream 
velocity  on  local  heat  transfer  rates. 

Yawed  cylinders  are  of  interest  to  aerodynamicis ts  because 
of  the  correlation  of  a  yawed  cylinder  to  a  swept  back  aircraft 
wing.   Researchers  and  designers  may  also  encounter  yawed  cylinders 
in  heat  exchangers  and  other  applications. 

No  information  is  currently  available  on  the  effect  of 
flow  oscillations  on  a  yawed  cylinder.   Viteri  [Ref.  10]   deter- 
mined the  average  heat  transfer  coefficients  of  yawed  cylinders 
in  uniform  flow,  and  Kraabel  [Ref.  11]   determined  the  local 
heat  transfer  coefficients  around  yawed  cylinders  in  uniform 
flow.   These  papers  provide  a  basis  for  validation  of  the  yawed 
cylinder  results  reported  here. 

It  has  been  demonstrated  by  Miller  [Ref.  2]       that  augmented 
convective  heat  transfer  can  be  attributed  to  the  reduction  of 
laminar  to  turbulent  transition  Reynolds  numbers.   Despard 
[Ref.  12]   showed  that   freestream  oscillations  promote  laminar 
boundary  layer  separation  in  the  presence  of  an  adverse  pressure 
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gradient.   Banning  [Ref.  13]   produced  data  that  revealed  that 
the  starting  vortex  associated  with  periodic  separation  from 
an  airfoil  in  oscillating  flow  produced  an  unexpectedly  strong 
boundary  layer  reattachment  which  markedly  improved  the  aero- 
dynamic performance  of  the  foil. 

In  summary,  the  effects  of  early  laminar  to  turbulent  boun- 
dary layer  transition  and  the  destabilized  boundary  layer  pro- 
duced slightly  augmented  heat  transfer  while  periodic  vortices 
shed  from  the  separation  point  can  cause  periodic  reattachment 
of  the  hydrodynamic  boundary  layer.   This  can  be  expected  to 
greatly  enhance  the  heat  transfer  in  the  afterbody  region. 

The  purpose  of  this  investigation  is  to  determine  the  effect 
of  large  flow  oscillations  on  the  local  heat  transfer  coeffi- 
cients around  both  normal  and  yawed  cylinders.   Additionally  it 
is  desired  to  develop  a  correlation  formula  for  the  average 
Nusselt  number  in  oscillating  flow  that  includes  the  effect 
of  yaw  angle,  Reynolds  number,  and  flow  frequency  and  amplitude. 
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II.   DESCRIPTION  OF  APPARATUS 

A.   MODEL 

The  model  used  in  the  investigation  was  a  four  inch  dia- 
meter cylinder.   It  was  fabricated  from  one  half  inch  wall 
thickness  fiberglass-phenolic  tubing  two  feet  long.   The  12 
inch  mid-span  section  was  instrumented  with  16  platinum  heater/ 
resistance  thermometers  oriented  parallel  to  the  axis  of  the 
cylinder  and  spaced  12  degrees  apart  from  center-to-center.   A 
sketch  of  the  model  showing  the  orientation  of  the  heater  strips 
is   shown  in  Figure  1. 

The  platinum  strips  used  were  12  inches  long  by  one-quarter 
of  an  inch  wide  by  one- thousandth  of  an  inch  thick.   Each  strip 
subtended  an  arc  of  7.16  degrees  of  the  circumference  of  the 
cylinder.   The  strips  were  spot  welded  to  brass  studs  that  were 
machined  flush  with  the  cylinder  on  the  outside  and  protruded 
into  the  interior  of  the  model  for  attachment  of  the  voltage 
sensing  and  current  carrying  wires.   Figure  2  shows  a  cross 
section  of  the  cylinder  wall  through  the  studs  and  strips. 
Voltage  drop  across  the  platinum  strip  was  measured  using  the 
voltage  sensing  wires  attached  to  each  stud.   Platinum  was 
used  to  both  heat  the  cylinder  and  determine  the  surface  tem- 
perature of  the  cylinder  because  of  its  large  linear  resistivity- 
temperature  relationship.   The  resistance  of  the  strips  was 
determined  from  the  voltage  drop  across  the  strip  which  was 
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then  to  be  used  to  find  the  temperature.   The  amount  of  heat 
convected  at  each  location  was  proportional  to  the  power  dis- 
sipated by  the  strip  which  was  equal  to  the  product  of  the  cur- 
rent squared  and  the  strip  resistance. 

The  entire  outside  of  the  cylinder,  including  the  platinum 
strips,  was  covered  with  a  thin  layer  of  hard  epoxy.   The  hard 
epoxy  layer  was  machined  and  polished  on  a  lathe  to  reduce  the 
thickness  over  the  strips  to  less  than  one-hundredth  of  an  inch 
and  to  produce  a  smooth  heat  transfer  surface.   A  photograph 
of  the  completed  model  showing  the  mounting  stud  and  bracket 
and  electrical  leads  is  shown  in  Figure  3. 

Two  copper-cons tantan  thermocouples  were  mounted  oji  the 
inside  surface  of  the  cylinder  at  90  degrees  and  270  degrees 
from  the  front  stagnation  point.   These  thermocouples  were  used 
to  monitor  the  interior  wall  temperature  of  the  cylinder  in  order 
to  estimate  the  heat  loss  to  the  inside  of  the  cylinder. 

B.   WIND  TUNNEL 

The  experiment  was  carried  out  in  the  oscillating  flow  wind 
tunnel  of  the  Naval  Postgraduate  School.   Photographs  of  the 
facility  are  shown  in  Figures  4  and  5.   The  oscillating  flow 
wind  tunnel  is  a  low  speed,  open  circuit  wind  tunnel  with  the 
inlet  bell  extending  out  of  the  building.   The  tunnel  inlet 
bell  is  eight  feet  square  and  the  test  section  is  two  feet 
square,  providing  a  16.1  contraction  ratio.   Three  high  solidity 
screens  located  in  the  inlet  section  upstream  of  the  nozzle 
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produce   freestream  turbulence  intensities  of  less  than  0.5 
percent  of  the   freestream  velocity  for  uniform  flow. 

The  wind  tunnel  drive  consists  of  two  axial  flow  fans  in 
series,  each  of  which  has  an  internal,  100  horsepower,  direct 
connected,  1750  RPM  motor.   The  fan  blades  are  internally  adjust- 
able through  a  pitch  range  of  25  to  55  degrees,  providing  a 
wide  operating  base.   A  set  of  variable  inlet  vanes,  located 
immediately  upstream  of  each  fan,  are  externally  adjusted  to 
provide  control  of  test  section  velocity.   These  vanes,  of  radial 
configuration,  preswirl  the  air  in  the  direction  of  fan  rotation 
to  reduce  fan  capacity.   The  total  range  of   freestream  velo- 
cities possible  is  from  ten  to  250  feet  per  second. 

A  plan  view  sketch  of  the  wind  tunnel  showing  the  location 
of  the  fans  and  inlet  vanes  in  relation  to  the  test  section  is 
shown  in  Figure  6. 

Two  fundamental  methods  are  available  to  superimpose  oscil- 
lating flow,  on  a  mean  velocity.   Both  Nickerson  [14]   and 
Hori  [15]   introduced  oscillations  by  oscillating  their  models 
in  a  steady  flow  environment.   This  method  severely  restricts 
the  range  of  attainable  frequencies  because  of  mechanical  com- 
plications and  instrumentation  difficulties. 

The  other  approach  is  to  provide  an  oscillating  flow  environ- 
ment over  a  stationary  model.   Hill  [16]   used  a  sliding  shutter 
to  impose  oscillations  on  the   freestream  but  was  limited  to 
low  frequencies  because  of  mechanical  difficulties.   Feiler 
and  Yeager  [17]   used  a  siren  mounted  upstream  of  an  eight  inch 
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diameter  wind  tunnel  to  expose  small  models  to  frequencies 
ranging  from  34  to  300  Hz.  with  root-mean-squared  amplitudes 
of  up  to  65  percent  of  the  mean  velocity.   The  most  successful 
method  of  obtaining  oscillating  flows  with  a  large  range  of 
frequencies,  amplitudes  and  flow  velocities  in  a  moderate  sized 
tunnel  is  that  employed  by  Karlson  [18]   and  later  by  Miller 
[19].   A  rotating  shutter  valve,   immediately  downstream  of  the 
test  section  is  used  to  superimpose  a  periodic  variation  of 
velocity  on  the  mean  flow.   The  method  used  in  this  investiga- 
tion is  identical  to  that  used  by  Miller. 

The  shutters  consist  of  four  horizontal,  equally  spaced, 
shafts  mounted  in  a  plane  just  downstream  of  the  test  section. 
The  shafts  are  slotted  to  receive  blades  of  various  widths, 
forming  a  set  of  four  butterfly  valves  spanning  the  wind  tunnel 
Figure  7  is  a  picture  of  the  shutter  valve  assembly.    The 
shutters  are  driven  synchronously  by  a  pulley  and  belt  arrange- 
ment and  a  vari-drive  motor  and  gearcase,  and  can  be  varied 
from  0.1  -  250  Hz.   The  amplitude  of  the  oscillation  is  a  func- 
tion of  the  width  of  shutter  blades,  frequency  of  rotation  of 
the  shutters  and  the   freestream  velocity.   Five  sets  of  blades 
are  available  producing  tunnel  blockages  of  33,  50,  66,  83 
and  100  percent  of  the  flow  area  in  the  fully  closed  position. 
In  this  investigation  blades  providing  66  and  100  percent 
blockage  were  used  resulting  in  amplitudes  of  5  to  69  percent 
of  the  mean   freestream  velocity. 
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The  test  section  of  the  tunnel  is  fabricated  from  two-two 
inch  thick,  25  inch  wide  by  18  feet  long  pieces  of  aluminum 
as  the  upper  and  lower  walls  and  three-two  inch  thick  stress 
relieved  lucite  panels  on  each  side  as  side  walls.   The  front 
panels  are  hinged  at  the  top  and  can  be  hydraulically  opened 
and  closed  for  access  to  the  test  section. 

For  this  investigation  the  rear  center,  lucite  panel  has 
been  replaced  by  a  two  inch  thick  plywood  panel  to  facilitate 
mounting  the  model.   The  model  was  mounted  by  a  chordwise 
extension  held  in  a  rotatable  clamp  mounted  flush  in  the  center 
of  the  plywood  panel.   The  other  end  of  the  model  is  supported 
by  a  3/8  inch  thick  "Y"  shaped,  faired  bracket,  that  is  affixed 
to  the  tunnel  ceiling  and  floor,  and  rests  against  the  front 
lucite  panel.   The  model  is  shown  in  place  with  the  "Y"  bracket 
clearly  visible  in  Figure  8.   Electrical  leads  were  taken  out 
through  the  plywood  back  panel  where  they  were  connected  to  the 
power  supplies  and  the  voltage  drop  sensing  circuits.   The 
model  was  yawed  by  affixing  extension  pieces  to  each  end  that 
were  angled  and  sized  to  completely  span  the  tunnel  at  each 
angle  of  yaw. 

Power  to  heat  the  strips  was  controlled  by  sixteen  variacs 
that  powered  16  individual  DC  power  supplies  (see  working  cir- 
cuit, Figure  9).   A  precision  shunt  (R=0.001  )  was  placed  in 
series  with  each  strip  and  the  voltage  drop  across  the  known 
resistance  was  used  to  determine  the  current.   The  voltage 
drop  across  each  strip  was  measured  using  the  potential  sensing 
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leads  that  were  attached  to  the  studs  on  the  inside  of  the  model 
Instantaneous  velocity  measurements  for  the  purpose  of  free- 
stream  turbulence  determination  in  uniform  flows  and  amplitude 
and  frequency  measurement  in  oscillating  flows  were  made  using 
a  linearized  constant  temperature  hot-wire  anemometer  described 
in  Reference  20.   An  oscilloscope  and  a  Ballantine  True  R.M.S. 
meter  were  used  to  monitor  hot-wire  signals.   Mean  velocity  of 
the  free  stream  was  measured  using  a  Prandtl  Probe  and  a  water- 
filled  micromanome ter . 

Figure  10  shows  a  pair  a  typical  velocity  wave  forms .   The 
vertical  axis  is  velocity  and  the  horizontal  axis  is  time.   The 
roughness  of  the  upper  tracing  is  typical  of  low  amplitude  oscil- 
lations which  are  significantly  more  difficult  to  produce  in 
true  sinusodial  form. 

The   freestream  temperature  and  the  model  internal  wall 
temperature  were  measured  using  copper-cons tantan  thermocouples 
and  a  Newport  digital  temperature  indicator. 

The  frequency  of  oscillation  was  also  monitored  with  an 
optical  chopper  driven  by  a  shutter  valve  shaft  and  a  digital 
counter . 
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III.   EXPERIMENTAL  PROCEDURE 

A.   CALIBRATION  OF  MODEL 

The  platinum  strips  were  calibrated  as  resistance  thermo- 
meters by  determining  their  electrical  resistance  as  a  function 
of  temperature  in  an  electric  oven.   Seven  separate  tempera- 
tures were  used,  evenly  spaced  between  ambient  (60  F)  and 
200°F.   Although  the  platinum  heater- thermome ter  strips  were 
extremely  pure  and  are  expected  to  have  a  linear  resistance- 
temperature  characteristic,  calibration  was  carried  out  to 
determine  the  exact  slope  and  intercept  of  each  strip. 

The  calibration  circuit  employed  (see  Figure  9)  was  powered 
by  a  low  impedance  DC  voltmeter  calibration  standard.   This 
proved  to  be  an  extremely  precise,  drift  free  power  supply.  The 
voltage  drop  across  the  precision  resistance  (1.0   +  0.15%)  was 
recorded  as  was  the  voltage  drop  across  the  heater  strip  itself 
The  current  was  then  computed  from  the  voltage  drop  across  the 
precision  resistance  divided  by  one  ohm,  and  the  resistance  of 
the  strip  for  that  temperature  was  computed  from  the  voltage 
drop  across  the  strip  divided   by  the  current.   Temperatures 
in  the  oven  were  monitored  with  a  Newport  Laboratories  Digital 
Temperature  Indicator  and  a  copper-cons tantan  thermocouple. 

After  determining  the  resistances  as  a  function  of  the 
temperature  a  straight  line  was  fitted  to  the  data  using  linear 
regression.   The  resulting  equations  were  used  throughout  the 
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subsequent  experimentation  and  data  reduction  to  determine  sur- 
face temperatures. 

B.   ADJUSTMENT  OF  SURFACE  HEATERS 

Before  each  series  of  runs  the  heater  strips  were  energized 
and  the  wind  tunnel  was  run  with  a  low  velocity  uniform  flow 
in  order  to  allow  the  model  substrate  to  come  to  thermal  equili- 
brium.  Once  this  was  accomplished  the  shutters  were  turned  on 
and  adjusted  to  the  desired  frequency  and  the  velocity  was 
adjusted  to  the  magnitude  corresponding  to  the  Reynolds  number 
desired.   At  this  point  the  model  was  again  allowed  to  come  to 
thermal  equilibrium. 

Once  this  initial  thermal  equilibrium  was  achieved  the 
voltage  drop  and  current  through  each  strip  was  recorded  and 
inserted  into  a  predictor-corrector  program    (see  Appendix  F)« 
This  program  calculated  a  new  required  voltage  drop  based  on 
desired  strip  temperature,  current  strip  temperature,  and  the 
known  resistance- temperature  calibration  curve  for  that  strip; 
the  strip  was  then  adjusted  to  the  new  desired  voltage  drop 
through  adjustment  of  its   power  supply.   Following  adjustment 
of  the  strips  the  model  was  again  allowed  to  return  to  thermal 
equilibrium. 

These  adjustments  were  performed  as  many  times  as  necessary 
until  all  the  strips  were  at  the  same  temperature,  ensuring  an 
isothermal  heat  transfer  surface  that  minimized  heat  conduction 
between  strips . 
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The  temperature  for  a  particular  run  was  determined  from 
a  sampling  of  the  strip  temperatures  before  the  first  set  of 
adjustments  was  made  and  was  selected  to  minimize  the  number 
of  adjustments  required.  In  every  case  however  the  temperature 
was  chosen  to  be  over  110  F  in  order  to  maintain  at  least  a 
50  temperature  difference  between  the  strip  and  the  ambient 
temperature . 

C.   FREESTREAM  CONDITIONS 

A  series  of  runs  was  made  with  both  the  45   yawed  model 
and  the  normal  model  covering  a  wide  range  of  operating  condi- 
tions.  Each  model  was  tested  with  both  the  six  inch  blades 
(100  percent  tunnel  blockage) and  the  four  inch  blades  (66  per- 
cent tunnel  blockage),  resulting  in  oscillation  amplitudes  from 
five  to  69  percent  of  freestream  mean  velocity.   Each  combina- 
tion of  frequency  and  amplitude  was  run  at  two  freestream  velo- 
cities corresponding  to  nominal  diametral  Reynolds  numbers  of 
300,000  and  150,000.   Frequencies  investigated  were  0,  5,  10, 
22,  50,  100  and  125  Hz. 


D.   DATA  REDUCTION 

The  data  reduction  technique  is  described  in  Appendix  A 
and  a  sample  calculation  is  presented  in  Appendix  B.   A  detailed 
energy  balance  is  necessary  in  order  to  determine  what  portion 
of  the  total  Joulian  heat  generated  by  platinum  resistance 
heater  is  transferred  from   the  surface  by  convection.   Part  of 
the  energy  was  conducted  through  the  thin  epoxy  coating  and 
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convected  directly  to  the  freestream  or  was  radiated  into  the 
environment.   The  balance  of  the  energy  was  conducted  into  the 
substrate  of  the  cylinder  where  some  was  conducted  and  con- 
vected into  the  center  of  the  cylinder,  some  was  lost  to  the 
ends  of  the  cylinder,  and  some  was  conducted  to  the  gaps  between 
the  strips  where  it  was  convected  to  the  freestream  or  was  radia- 
ted to  the  environment. 

The  energy  lost  from  the  substrate  side  of  the  surface 
heaters  is  hard  to  determine  analytically  since  it  consists  of 
both  conduction  to  the  intraheater  gaps  and  subsequent  loss  by 
convection  and  radiation  as  well  as  conduction  to  the  inner 
surface  of  the  model.   Therefore  a  two-dimensional  Teledeltos 
paper  model  was  constructed  to  determine  the  relationship 
between  the  thermal  resistances  to  the  inner  and  outer  surfaces. 
This  procedure  and  the  results  are  discussed  in  Appendix  D. 
The  results  of  the  Teledeltos  paper  experiment  were  placed  in 
equation  form  and  solved  interatively  for  the  external  covec- 
tion  heat  transfer  coefficient,  h,  for  each  strip. 

The  other  major  loss  of  energy  from  the  system  is  by  radia- 
tion to  the  environment.   These  losses  can  be  appreciable  and 
an  accounting  for  this  loss  was  developed  as  part  of  the  data 
reduction  program.   Details  are  described  in  Appendix  C. 
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IV.   RESULTS  AND  DISCUSSION 

A.   STEADY  FLOW 

Six  steady  flow  runs  were  made,  two  with  the  model  normal 

to  the  flow  at  Nn   =  150,000  and  300,000  and  four  runs  with 

Re 

the  model  at  45°  to  the  flow  at  NRe  =  35,000,  100,000,  150,000, 

and  300,000.   The  results  of  these  runs  are  shown  graphically 

in  Figures  13  through  18.   Results  of  the  model  at  90   (normal 

to  the  flow)  at  N_   =  150,000  are  compared  in  Figure  19  with 

Re 

the  theory  of  Froessling  [Ref.  24]   and  the  experimental 

results  of  Schmidt  [Ref.  25  1   for  N_   =  170,000.   The  discrepancy 

Re 

between  theory  and  both  of  the  experimental  results  is  expected 

and  is  a  consequence  of  the  effect  of  freestream  turbulence 

near  the  stagnation  point  [Ref.  26]. 

Comparison  of  the  results  of  the  present  experiment  with 

those  of  Schmidt  indicates  two  interesting  points.   The  two 

agree  within  20  percent  up  to  an  angle  of  150   from  the 

stagnation  point,  beyond  which  they  diverge  markedly,  differing 

by  70  percent  at  180  .   This  sharp  increase  in  heat  transfer 

above  that  noted  by  Schmidt  is  typical  of  all  the  results 

observed  in  the  present  investigation  and  shows  no  correlation 

with  frequency  number  or  Reynolds  number.   Agreement  within 

20  percent  is  good  however  and  this  lends  a  high  degree  of 

confidence  in  the  experimental  techniques  and  results. 
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No  local  heat  transfer  data  is  available  for  the  steady 
flow  case  with  45   yaw.   However,  the  data  of  Viteri  [Ref.  10] 
can  be  interpolated  to  arrive  at  an  average  Nusselt  number  of 
128.3  at  a  Reynolds  number  of  31,900.   At  a  Reynolds  number  of 
35,630  the  present  experiment  produced  an  average  Nusselt  number 
of  122.1.   This  close  agreement  lends  confidence  in  at  least 
the  average  Nusselt  Numbers  obtained  in  the  45   yawed  case. 


B.   NONSTEADY  FLOW 

Two  thirds  of  the  nonsteady  freestream  conditions  did  not 
produce  substantial  alteration  in  heat  transfer  when  compared 
to  the  steady  flow  case.   Figures  20  through  25  depict  the 

cases  with  the  90   model  that  show  little  or  no  altered  heat 

o 
transfer,  and  Figures  26  through  39  depict  those  of  the  45 

model.   In  each  of  these  figures  the  steady  flow  data  has  been 
plotted  along  with  the  nonsteady  data  for  comparison. 

Twelve  of  the  38  freestream  conditions  investigated  pro- 
duced substantially  altered  local  heat  transfer  rates.   These 
results  are  shown  graphically  in  Figures  40  through  49  for  the 
90   model,  and  Figures  50  and  51  for  the  45   model.   The 
changes  in  local  Nusselt  number  took  place  almost  exclusively 
in  the  afterbody  region  of  the  model,  and  resulted  from  altered 
separation  and  flow  patterns  in  the  wake  of  the  cylinder.   Fig- 
ure 40  is  a  typical  example  of  these  results  and  demonstrates 
enhanced  heat  transfer  in  the  wake  area  of  the  90   model.   The 
data  reported  in  Figures  48  through  51  exhibit  diminished  heat 
transfer  rates.   Figure  50  is  typical  of  this  class  of  data. 
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In  order  to  simplify  the  analysis  of  the  results,  the 
average  Nusselt  number  was  calculated  for  each  freestream  con- 
dition (see  Table  I).   Figures  52  and  53  report  average  Nusselt 
numbers  versus  Reynolds  number  for  the  90   model  and  the  45 
model  respectively.   In  the  90   case  (Figure  52)  a  curve  has 


been  drawn  representing  Nu  =  N 


0.5716 


Re 


and  the  Amplitude  and 


Storuhal  numbers  have  been  annotated  for  each  data  point.   Of 
interest  is  the  fact  that  data  with  higher  Amplitude  numbers 
are  those  which  tend  to  exhibit  augmented  average  heat  transfer 
rates.   Moreover,  higher  Amplitude  numbers  seem  to  correspond 
to  greater  heat  transfer  augmentation.   Although  the  data  is 
sparse,  a  pattern  of  larger  Nusselt  numbers  for  larger  Ampli- 
tude numbers  at  the  same  Reynolds  number  seems  to  exist.   Three 
apparent  anomalies  exist  associated  with  the  data  from  run  num- 
bers 15,  17  and  20.   It  is  possible  that  an  error  in  recording 
the  amplitude  data  was  made  and  that  the  Amplitude  numbers  for 
runs  17  and  20  were  reversed.   If  this  were  true  then  both 
runs  17  and  20  would  conform  to  the  apparent  pattern.   When  the 
Frequency  numbers,  Nusselt  numbers  and  Reynolds  numbers  of  runs 
15  and  17  are  compared,  it  is  seen  that  they  are  very  similar 
but  that  the  Amplitude  numbers  are  considerably  different.   If 
the  previously  noted  pattern  for  higher  Nusselt  number  for 
higher  Amplitude  number  is  true  then  the  Amplitude  number  of 
run  15  is  suspect  and  may  be  too  low. 


27 


Average  Nusselt  number  data  taken  for  the  45   model  is 
shown  in  Figure  53.   No  apparent  pattern  of  heat  transfer  altera- 
tion is  evident  in  this  plot  even  though  Amplitude  numbers  as 
high  as  0.44  and  0.65  appear  in  the  data. 

C.   NUSSELT  NUMBER  CORRELATION 

The  data  obtained  in  this  investigation  is  too  sparse  to 
allow  a  correlation  between  heat  transfer  rate  and  freestream 
conditions  to  emerge.   Only  a  possible  pattern  is  suggested  by 
the  results  and  more  data  covering  a  wide  range  of  Reynolds 
numbers  and  Amplitude  numbers  is  necessary  to  produce  such  a 
correlation . 
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V.   CONCLUSIONS 

Although  there  is  insufficient  data  to  draw  any  definite 
conclusions  or  to  derive  a  correlation,  several  trends  may  be 
observed.   -In  the  90   case,  significant  increases  in  the  local 
Nusselt  numbers  are  seen  in  the  wake  region  and  are  due  to 
changes  in  the  separation  and  wake  structure.   Moreover,  aug- 
mentation of  local  Nusselt  numbers  seems  to  increase  with 
increasing  frequency  as  suggested  by  the  data  shown  in  Figures 
40  through  47.   Increases  in  average  Nusselt  numbers  is  appar- 
ently related  to  the  Amplitude  numbers  as  discussed  above  and 
shown  in  Figure  52.   With  the  exception  of  the  three  points 
associated  with  runs  15,  17  and  20,  discussed  in  Chapter  IV, 
the  pattern  exhibited  is  for  higher  average  Nusselt  numbers  to 
be  associated  with  greater  amplitudes. 

In  the  case  of  the  45   yawed  cylinder,  no  general  patterns 
are  discernable  and  it  appears  that  local  Nusselt  numbers 
exhibit  no  significant  change  over  a  wide  range  of  frequencies 
and  amplitudes. 

This  investigation  has  served  to  initiate  a  study  of  heat 
transfer  from  a  cylinder  in  oscillating  flow.   Additionally 
it  has  validated  the  techniques  to  be  used  for  further  study. 
More  extensive  data  collection  will  be  required  before  defini- 
tive conclusions  can  be  drawn  or  correlations  made. 
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APPENDIX  A 


DATA  REDUCTION 


The  data  was  reduced  to  the  form  shown  in  Table  2  by  a 

specially  written  Fortran  computer  program  on  an  IBM  360  com- 
puter at  the  Naval  Postgraduate  School. 

The  raw  data  consisted  of  pitot-static  tube  4P  (inches 

HO),  freestream  temperature,  T  ,  (  F) ,  internal  model  tempera- 
Z  r 

tures  (heater  side  and  "other"  side,   F)  ,  frequency  of  oscil- 
lation (Hz),  amplitude  ratio  (%)  and  the  voltage  drop  across 
(volts)  and  current  passing  through  (amperes)  each  strip.   The 
total  power  dissipated  by  each  strip  is  simply: 

P  =  3.4149  x  E  x  I 

where  3.4149  is  a  conversion  factor  from  watts  to  BTU/Hr. 

From  the  measured  current  (I)  and  voltage  drop  (E)  each 
individual  strip  temperature  is  determined  by  comparing  its 
resistance  to  its  known  temperature  versus  resistance  charac- 
teristics.  The  resistance  of  a  strip  is: 


*-f 


From  the  linear  temperature-resistance  relationship  obtained 
during  calibration,  the  strip  temperature  can  be  obtained: 

R-C 


30 


where  C  and  S  are  the  known  ordinate  intercept  and  slope 
obtained  from  the  individual  resistance  versus  temperature 
calibration  curve. 

The  other  unknown  temperature  required  by  the  data  reduc- 
tion program  is  the  inner  wall  temperature,  T  .   The  inner  wall 
temperatures  are  known  at  two  points.   Thermocouples  are  mounted 
on  the  wall  at  locations  90   and  270   from  the  stagnation  point 
at  mid-span  (see  Figure  2).   The  inner  wall  temperature  is 
assumed  to  vary  linearly  between  these  two  measured  values. 

With  this  information  h,  the  heat  transfer  coefficient  may 
be  calculated.   The  following  formula  (developed  in  Appendix 
D)  is  used : 


s     F 


1  + 


D 


R 


T   -  T 
^_S _1 

R~ 


- 

RI 

h  is  contained  in  the  expressions  for  R  ,  R  /Rn,  and 
RT/R_  as  is  shown  in  Appendix  D.   This  equation  was  solved 
iteratively  to  obtain  h. 

With  the  local  h  for  each  strip,  the  local  Nusselt  number 
is  calculated  from: 


Nu  = 


hD 
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APPENDIX  B 


SAMPLE  CALCULATION 


For  the  purpose  of  this  sample  calculation  strip  Nr  6  from 
data  run  Nr  20  was  selected. 


The  information  recorded  is: 

Voltage  drop 

Current 

Pitot-Static  Pressure 

Ambient  Temperature 

Inside  Temperature: 
Heater  Side 
Other  Side 

Frequency 

Amplitude  Ratio 

Blades 

Model 

Barometric  Pressure 

Freestream  Mean  Velocity 


V  = 


2gHW  (2rwater)  -  1) 
&    v  yair  '  ' 


] 


1/2 


=  134.4 


ft 
sec 


E  =  1.62  Volts 
I  =  7.47  Amps 


HW  =  4.16  inches  Yi    Q 


TA  =  64.0  F 


TIHTR  =  105. 0°F 
TIOTH  =   77.0°F 


Freq  =  50.0  CPS 

N.  =  0.31 

A 

Medium  (4") 

_  o 

90   to  stream 

30.04  inches  Hg 


62.4 


2(32.2  %(4.16  in  H,0)  (-A  fe 

_ Z  ILL.     xn 


lbf 
ft3 


0.07654 


lbf 
ft3 


-  1 


1/2 
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Frees tr earn  Reynolds  number 


ft 


Re    V 


134.8  —  (0.333  ft) 

s 

0.000166  ft2/s 


=  269,800 


Temperature  of  the  strip 


T  = 


!-c 


E  m    1.622  Volts 
I    7.47  Amps 


=  0.2171  ohms 


C  =  0. 173955807 


S  =  0.000390625 


Therefore , 


T  = 


0.  2171-0. 173955807 
0.000390625 


=  110.5  F 


Total  Heat  Dissipated  by  the  strip: 

qs  =  EI  =  (3.4149  —-jjB----)  (1.622  Volts)(7.47  Amps) 


-  41.38  Btu 


Also  : 


D  I 


(Ts  -  TF)(1  +I£)  +  (Ts 


v^"1 


R, 


k  +  ih. 
khTA 


h 

R, 


al  +  a2  hT 
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—  =  b1  hT  +  b2  hT    +  b3  hT 


After  substitution,  the  equation  may  be  solved  for  h  . 
(the  constants  a..  ,  a„,  b,,  b?    and  b ,,  are  given  in  Appendix  D) 
By  iteration: 


h   =  36.9 


Btu 


hr  -  ft2  -  °F 


The  radiant  heat  transfer  coefficient 


hr  =  e<T(Ts  +  iA)(is2  +  ta2) 


Assuming  T   is  the  same  as  the  fluid  temperature: 

h   =  (0.73) (0.1718  x  10~8)(570.5  +  524.0) 


x  (570. 52  +  524. 02) 


=  0.824 


Btu 


hr  -  ft2  -  °F 


Therefore,  the  local  convective  heat  transfer  coefficient  is 


h  =  b   -  h 
T     r 


=  36.9  -  0.824 
Btu 


=  36.1 


hr  -  ft2  -  °F 
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The  local  Nusselt  number 


Nu  = 


hd  =  (36.1)  (0.333) 
k      0.014738 


=  815.9 
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APPENDIX  C 


RADIATION 


For  a  uniformly  convex  gray  body  with  emissivi ty  £■  ,  at  a 
temperature  T,  radiating  into  an  enclosure  of  temperature  T— , 
the  net  radiant  heat  exchange  is: 


qr  =   <f€  A(T4  -  t/) 


where  <5"*  is  the  S  tef  an-Bo  ltzmann  constant,  and  A  is  the  surface 
area  of  the  emitter. 

The  above  equation  can  be  put  in  the  same  form  as  the  con- 
vection equation: 


q   =  h  A(T  -  T  ) 


where  h   is  the   effective  radiation  heat  transfer  coefficient 

r 


given  by: 


hr  =  6r«-   (T  +  Tjd2  +  T^2) 


Taking  £  =  0.73  for  epoxy  coating  over  aluminum  [Ref.  22] 
and  the  following  typical  temperatures: 
T  =  120°F  =  580°R 
Tm=    60°F   =  520°R 
Gives : 
h   =  0.73  (0.1714) (10-8)  (580  +  520)(5802  +  5202) 


=  0.835 


Btu 

hr -ft2-F 
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This  value  falls  within  two  to  10  percent  of  the  total 
heat  transfer  coefficient  and  therefore  cannot  be  neglected. 
Therefore,  the  convective  heat  transfer  coefficient,  h,  is 
obtained  by  subtracting  the  radiant  heat  transfer  coefficient, 
h  ,  from  the  experimentally  derived  total  heat  transfer  coef- 
ficient ,  that  is , 


h  -  hT  -  hr 
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APPENDIX  D 
CONDUCTION  LOSSES 

Although  the  thermal  resistance  of  the  fiberglass  substrate 
is  very  high  compared  to  the  thermal  resistance  between  the 
strip  and  the  fluid,  some  heat  is  lost  to  the  interior  of  the 
model.   This  heat  loss  must  be  deducted  from  the  total  heat 
input  in  order  to  arrive  at  the  actual  heat  transferred  to  the 
air.   Heat  loss  due  to  conduction  in  the  longitudinal  axis  direc- 
tion and  radially  through  the  small  brass  studs  was  neglected. 

A  Teledeltos  paper  model  was  used  to  estimate  the  heat  loss 
to  the  interior  of  the  model.   Symmetry  was  used  to  reduce  the 
size  of  the  required  model  (see  Figure  11). 

Teledeltos  paper  is  a  graphite  covered  paper  having  a  nomi- 
nal resistivity  of  1000  ohms  per  inch  square  and  is  highly 
isotropic.   Teledeltos  paper  is  extremely  useful  for  two- 
dimensional  heat  conduction  modeling  in  cases  such  as  this 
which  have  boundary  conditions  and  a  physical  shape  that  are 
inconvenient  for  an  exact  solution. 

The  heater  strip  itself  is  simulated  with  a  strip  of  silver 
conductive  paint  that  is  placed  on  the  paper  in  exactly  the 
same  position  as  the  platinum  strip  is  on  the  actual  model. 
The  inside  of  the  model  segment  is  considered  isothermal  so 
the  analogous   surface  of  the  Teledeltos  paper  model  was 
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held  at  uniform  potential  with  a  conductive  strip  of  silver 
paint.   The  left  and  right  hand  adiabatic  boundaries,  between 
strips,  are  easily  modeled  by  leaving  them  bare.   The  convective 
boundary  at  the  gap  between  strips  is  modeled  by  allowing  the 
Teledeltos  paper  to  extend  beyond  the  actual  cylinder  surface 
for  a  distance  to  simulate  the  convection  resistance.   The 
required  distance  is  determined  by: 


L  -E  s 


where  k  is  the  conductivity  of  the  actual  substrate,  h,  is  the 
convective  coefficient  being  modeled  and  S  is  the  scale  factor, 
in  this  case,  S  =  48.   Additionally  the  excess  length  was  slit 
into  eight  equal  width  strips  lengthwise  to  force  the  flux 
lines  parallel  to  each  other  and  perpendicular  to  the  surface. 
In  order  to  model  h's  of  10,  20,  50  and  100  Btu/hr • f t2 • F ,  L's 
of  14.4,  7.2,  2.88  and  1.44  inches  respectively  were  used. 

A  strip  temperature  of  120  F  with  an  inner  temperature 
of  100  F  and  a  fluid  temperature  of  60  F  were  used.   These 
temperatures  were  scaled  to  appropriate  voltages  that  would 
avoid  excessive  Joulian  heating  of  the  Teledeltos  paper. 

In  operation  the  voltages  were  measured  at  17  equally 
spaced  intervals  along  the  line  corresponding  to  the  cylinder 
surface  and  at  five  equally  spaced  points  along  a  line  cor- 
responding to  0.08  inches  from  the  inner  surface.   This  allows 
the  determination  of  the  ratio  of  heat  loss  to  the  inside  of 
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the  model,  compared  to  the  heat  loss  indirectly  through  the 
air  (via  the  substrate). 

Considering  the  lumped  electrical  resistance  network 
analogous   to  the  heat  transfer  problem  (see  Figure  12)  it  is 
noted  that  the  power  input  to  the  strip  is  equal  to  the  sum  of 
the  losses  to  the  air,  to  the  inside  of  the  cylinder  to  the 
gap  between  the  strips,  or: 


qD  +  q0  +  qI 


since  in  general: 


q  = 


4  T 


R 


TH 


this  becomes : 


TS  -  TF    TS  -  TF    TS  -  TI 
q   =  4.  +  

qs      RD  RQ  R 


Where  T  ,  T   and  T,  are  the  temperatures  of  the  strip,  the 

fluid  and  the  inside  respectively  and  R  ,  R   and  RT  are  the 

thermal  resistances  from  the  strip  to  the  air,  from  the  strip 

through  the  gap  to  the  air  (R^  =  R   +  R,  )  and  through  the 

u     a     b 

substrate  to  the  inside  of  the  cylinder  respectively  (see 
Figure  12 )  . 

This  equation  can  be  manipulated  into  the  following  form: 


T   -  T 
S     F 


1  + 


D 


+ 


T   -  T 
S      I 
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The  resistance  to  heat  loss  from  the  strip  directly  to  the 
air,  R  ,  is  due  to  the  thin  layer  of  epoxy  on  the  strip  and 
the  convection  film  resistance.   The  conduction  resistance 
due  to  the  epoxy  is  equal  to  -£/kA  where  JL   is  the  thickness  of 
the  epoxy  layer,  k  is  the  thermal  conductivity  of  the  epoxy 
and  A  is  the  surface  area  of  the  strip.   The  convective  resis- 
tance is  equal  to  1/hA,  where  h  is  the  convection  heat  transfer 
coefficient.   Summing  these  we  obtain: 

8   =  JE  +  -^h 

D     khA 

The  quantities  R  /R   and  R  /R   are  derived  experimentally 
from  the  Teledeltos  paper  analogue  and  are  in  the  form  of 
polynomial  equations  in  h: 


R0 

■~   =  1.546744  +  0.0320672  h 

RD 

RI  2  3 

=■=  =  0.106105h  +  0.0047559h   +  0.0000079h 

RD 

These  resistances  were  used  in  the  heat  transfer  equation 
to  solve  for  h  (see  Appendix  B) . 
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APPENDIX  E 
UNCERTAINTY  ANALYSIS 

For  the  purpose  of  this  analysis,  the  method  of  Kline  and 
McClintock.  [Ref.  23]   will  be  used,  assuming  that  all  variables 
have  a  normal  Gaussian  distribution.   Possible  dimensional  error 
in  strip  size,  placement  and  other  dimensions  were  considered 
so  small  as  to  be  negligible.   The  tabulated  values  of  k  and 
e  were  assumed  to  be  correct.   All  temperatures  measured  with 
thermocouples  were  assumed  to  have  an  uncertainty,  due  to 
readout  error,  of  0.1  degree  F. 

The  strip  temperature  uncertainties  derive  from  the 
+0.003  mV  uncertainty  in  reading  the  digital  voltmeter  used  in 
calibration   and  the  +0.002  V  and  +0.01  amp  uncertainty  in 
reading  the  meters  during  data  runs. 

The  uncertainty  in  measured  resistances  during  calibration 
is 


A  R  =  R 


2         2 


1  1/2 


=  +0.00031  ohms 

The  resistances  obtained  in  calibration  are  used  to  obtain 
a  straight  line  equation  for  temperature  vs  resistance. 

R  =  S-T  +  C 
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The  uncertainty  of  C  is  assumed  to  be  of  the  same  order  as 
R.   The  uncertainty  of  S  is: 


A  S  =  S 


(±f)2  +  (^f)2 


-i  1/2 


=  +0.00000064 


Since  the  strip  temperature  during  a  data  run  is  given  by 


T  = 


f-« 


The  uncertainty  of  T  is  given  by 


AT       =    T 
S  S 


»    2  ,  T    2  .„    2 


-    1/2 


+    T 


T  L^}2  +   (z|s}2 


-I  1/2 


=    +0.54WF 


The  values  of  R  /R   were  obtained  from  a  Teledeltos  paper 
model.   Uncertainties  in  measuring  a  potential  difference  between 
any  two  points  on  the  Teledeltos  paper  was  +0.02  volts.   This 
corresponds  to  a  uncertainty  in  temperature  of  +0-2  .   For 
RQ/RD  is  given  by: 


R 


ADT 


R     SADT 
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where  ADT  is  the  direct  convection  surface  area  times  (T.  -  T„) 

S     F 

and  SADT  is  the  indirect  convection  area  times  (T  local  -  T  ). 

r 

(T  local  is  determined  from  the  Teledeltos  paper  model.) 
The  uncertainty  of  SADT  is  given  by: 

SADT  =  SADT  (A^T)) 

=  +0.0011 


Similarly  for  ADT 


A  ADT    =    ADT     (^|) 


=    +0.0021 


Now 


R 


D 


(AADT.2  ASADT   '2 


"I    1/2 


ADT 


SADT 


=    +0.0088 


R_/R       is    given    by 


R. 


(T      -    T   ) 

ADT       h   US  V 


IADT 


k  (^    -    TF) 


Ay 


The  uncertainty  of  IADT  is  given  by 


A IADT  = 


IADT  (^f^-) 


=  +0.0  32 
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Thus  the  uncertainty  of  R  /R   is  given  by: 


R^    R_ 
R~    R 


D     D  L 


(M2I}2  +  (4IMI)2  +  (^I__)2  +  (^l   )2 


ADT 


IADT 


T  -T 
S   I 


T  -T 
S   F 


1/2 


=  +0.0  24 


q  is  given  by : 


a    =  EI 


The  uncertainty  in  q  is  given  by: 


Aq    =    q 


1/2 


(^f)2    +    i^)2 


=    +0.052 


Btu 
hr 


Because  of  the  complexity  of  the  formula  used  to  determine  h, 
h  was  found  by  iteration  of  the  following  implicit  equation: 


[ 


*  =^hil(Ts     V'1  +  aft 


)  +  (Tg  -  Tz) ( 


Blh  +  B2h2  +  B3h3 


The  uncertainty  of  h  was  determined  by 


Ah    = 


2  2  2 

AT  AT  AT 


2T, 


2T. 


2T 


1    1/2 


K     (Rn}     2(RjRn)}  ^     CRJ     2(RT/R^); 


V     2(R0/RD 
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The  partial  derivatives  were  evaluated  numerically  and  were 
found  to  be : 


2h 


2T 


=  -0.9137 


2h 


2T. 


=  1. 375 


2h 


2T. 


=  0.0698 


2h 


2(RQ/RD) 


-  4.  7619 


2h 


2(R1/RJ)) 


=  0.01 


therefore 


Ah    =  +0.51 


Btu 


Hr • ft  -F 


The  mean  velocity  is  given  by: 


V   =  [2gH  (g  Water  -1) 
M  a       air 


1   1/2 


The  only  parameter  here  that  is  open  to  uncertainty  (due  to 
the  assumption  that  property  values  are  correct)  is  the  H, 
or  height  of  water  column  in  the  micromanome ter .   Since 
AH  =  +0. 02  inches  : 

V, 


AV 


M 


M     ,4H. 
2     K    H; 


+0.324    ft/s 
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since : 


Re 


VD 
V 


A  N 


Re 


V 


=  +635. 


Nusselt  number  is  by  definition 


Nu  = 


hD 


therefore : 


A  Nu  =  Nu  (~) 

n 


=  +11.3 
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APPENDIX  F 
PREDICTOR  -  CORRECTOR  EQUATION 

In  order  to  simplify  and  shorten  the  process  of  adjusting 
strip  voltages  to  obtain  an  isothermal  surface  as  described 
in  Chapter  3,  a  predictor-corrector  program  was  developed  and 
implemented  on  a  TI-59  with  PC-100C  thermal  printer.   The  pro- 
gram listing  is  included  at  the  end  of  this  Appendix. 

The  program  requires  entry  of  the  ambient  temperature  and 
desired  tempera ture of  the  strips.   Then  the  strip  number  is 
entered  and  the  program  prompts  for  initial  voltage  drop  and 
initial  current  of  the  strip.   The  program  computes  strip 
temperature  using  the  formula: 


T  = 


!-' 


where  C  and  S  are  taken  from  the  calibration  curves. 

Next  the  program  calculates  the  voltage  drop  that  would  be 
expected  at  the  desired  strip  temperature  using: 

1  1/2 


'D 


R  Vi'W 


W  + 


Ix 
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Where  R  =  required  resistance,  from  calibration 

T  =  Fluid  temperature 

r 

S  =  Slope  of  calibration  curve 


The  variac  on  the  test  rig  is  then  adjusted  to  produce 
this  voltage  drop  across  the  strip.   The  strip  must  be  allowed 
to  come  to  thermal  equilibrium  after  an  adjustment  is  made. 
The  calibration  curves  assume  thermal  equilibrium  and  are  mean- 
ingless when  temperature  transients  are  taking  place. 

The  program  listed  requires  204  program  storage  locations 
and  76  memory  storage  registers.   The  TI-59  must  be  reparti- 
tioned  to  obtain  adequate  storage  registers  by  entering  8  (2nd) 
OP  17.   This  provides  319  program  storage  locations  and  79 
data  registers.   All  of  the  data  listed  must  be  entered  to  take 
advantage  of  the  full  capability  of  the  program. 

After  the  program  and  stored  data  are  loaded  the  ambient 
temperature  is  specified  by  entering  the  temperature  and  pres- 
sing key  "B".   The  desired  strip  temperature  can  then  be 
entered  and  the  key  "C"  depressed.   The  program  is  now  ready 
for  use . 

In  use  the  number  of  the  desired  strip  is  entered  and  key 
"A"  depressed.   The  printer  will  then  prompt  "E?"  and  the  initial 
voltage  drop  should  be  entered  and  the  "R/S"  key  depressed.   The 
program  will  echo  the  current  on  the  printer  and  will  proceed 
to  print  the  current  temperature  and  desired  voltage  drop.   The 
program  will  then  pause  and  is  ready  for  the   next  strip  number. 


49 


FIGURE  1  -  SKETCH  OF  THE  MODEL 
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FIGURE   2   -   CROSS-SECTIONAL  SKETCH  OF  THE  MODEL 
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CALIBRATION  CIRCUIT 
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FIGURE  9  -  SKETCH  OF  CALIBRATION  AND  WORKING 
CIRCUITS 
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FIGURE  10  -  TYPICAL  VELOCITY  WAVEFORMS 
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FIGURL  11  -  SKETCH  OF  TELEDELTOS  PAPER  MODEL 
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FIGURE  12  -  SKETCH  OF  THE  LUMPED  RESISTANCE 
MODEL 
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FIGURE  21  -  LOCAL  HEAT  TRANSFER  COEFFICIENTS  IN  OSCILLATING 
FLOW  THAT  SHOWED  UNALTERED  RESULTS 
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FIGURE  22  -  LOCAL  HEAT  TRANSFER  COEFFICIENTS  IN  OSCILLATING 
FLOW  THAT  SHOWED  UNALTERED  RESULTS 
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FIGURE  23  -  LOCAL  HEAT  TRANSFER  COEFFICIENTS  IN  OSCILLATING 
FLOW  THAT  SHOWED  UNALTERED  RESULTS 
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FIGURE  24  -  LOCAL  HEAT  TRANSFER  COEFFICIENTS  IN  OSCILLATING 
FLOW  THAT  SHOWED  UNALTERED  RESULTS 
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FIGURE  25  -  LOCAL  HEAT  TRANSFER  COEFFICIENTS  IN  OSCILLATING 
FLOW  THAT  SHOWED  UNALTERED  RESULTS 
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FIGURE  26  -  LOCAL  HEAT  TRANSFER  COEFFICIENTS  IN  OSCILLATING 
FLOW  THAT  SHOWED  UNALTERED  RESULTS 
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FIGURE  27  -  LOCAL  HEAT  TRANSFER  COEFFICIENTS  IN  OSCILLATING 
FLOW  THAT  SHOWED  UNALTERED  RESULTS 
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FIGURE  28  -  LOCAL  HEAT  TRANSFER  COEFFICIENTS  IN  OSCILLATING 
FLOW  THAT  SHOWED  UNALTERED  RESULTS 
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FIGURE  29  -  LOCAL  HEAT  TRANSFER  COEFFICIENTS  IN  OSCILLATING 
FLOW  THAT  SHOWED  UNALTERED  RESULTS 
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FIGURE  30  -  LOCAL  HEAT  TRANSFER  COEFFICIENTS  IN  OSCILLATING 
FLOW  THAT  SHOWED  UNALTERED  RESULTS 
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FIGURE  31  -  LOCAL  HEAT  TRANSFER  COEFFICIENTS  IN  OSCILLATING 
FLOW  THAT  SHOWED  UNALTERED  RESULTS 
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FIGURE  32  -  LOCAL  HEAT  TRANSFER  COEFFICIENTS  IN  OSCILLATING 
FLOW  THAT  SHOWED  UNALTERED  RESULTS 
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FIGURE  33  -  LOCAL  HEAT  TRANSFER  COEFFICIENTS  IN  OSCILLATING 
FLOW  THAT  SHOWED  UNALTERED  RESULTS 
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FIGURE  34  -  LOCAL  HEAT  TRANSFER  COEFFICIENTS  IN  OSCILLATING 
FLOW  THAT  SHOWED  UNALTERED  RESULTS 


83 


o 
oo 

C\J~ 


CD 
C\l' 


UJ 

az  o 


o 
r- 


o 
o 


O    -    UNIFORM    FLOW 
*    -    22.0     CPS 


* 


m      *      * 


6 


°0.00  36.00  72.00  108.00  144.00  180.0 


THETR     (DEGREES) 


45    DEGREE    MODEL 
RUN    NR    21 


REYNOLDS  NR   =155739 

FREQUENCY  NR  =  3.866  x  10-6 

AMPLITUDE  NR  =  0.040 

FIGURE  35  -  LOCAL  HEAT  TRANSFER  COEFFICIENTS  IN  OSCILLATING 
FLOW  THAT  SHOWED  UNALTERED  RESULTS 
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FIGURE  36  -  LOCAL  HEAT  TRANSFER  COEFFICIENTS  IN  OSCILLATING 
FLOW  THAT  SHOWED  UNALTERED  RESULTS 
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FIGURE  37  -  LOCAL  HEAT  TRANSFER  COEFFICIENTS  IN  OSCILLATING 
FLOW  THAT  SHOWED  UNALTERED  RESULTS 
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FIGURE  38  -  LOCAL  HEAT  TRANSFER  COEFFICIENTS  IN  OSCILLATING 
FLOW  THAT  SHOWED  UNALTERED  RESULTS 
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FIGURE  39  -  LOCAL  HEAT  TRANSFER  COEFFICIENTS  IN  OSCILLATING 
FLOW  THAT  SHOWED  UNALTERED  RESULTS 
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FIGURE  40  -  LOCAL  HEAT  TRANSFER  COEFFICIENTS  IN  OSCILLATING 
FLOW  THAT  SHOWED  ENHANCED  HEAT  TRANSFER 
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FIGURE  41  -  LOCAL  HEAT  TRANSFER  COEFFICIENTS  IN  OSCILLATING 
FLOW  THAT  SHOWED  ENHANCED  HEAT  TRANSFER 
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FIGURE  42  -  LOCAL  HEAT  TRANSFER  COEFFICIENTS  IN  OSCILLATING 
FLOW  THAT  SHOWED  ENHANCED  HEAT  TRANSFER 
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FIGURE  43  -  LOCAL  HEAT  TRANSFER  COEFFICIENTS  IN  OSCILLATING 
FLOW  THAT  SHOWED  ENHANCED  HEAT  TRANSFER 
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FIGURE  44  -  LOCAL  HEAT  TRANSFER  COEFFICIENTS  IN  OSCILLATING 
FLOW  THAT  SHOWED  ENHANCED  HEAT  TRANSFER 
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FIGURE  45  -  LOCAL  HEAT  TRANSFER  COEFFICIENTS  IN  OSCILLATING 
FLOW  THAT  SHOWED  ENHANCED  HEAT  TRANSFER 
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FIGURE  46  -  LOCAL  HEAT  TRANSFER  COEFFICIENTS  IN  OSCILLATING 
FLOW  THAT  SHOWED  ENHANCED  HEAT  TRANSFER 
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FIGURE  47  -  LOCAL  HEAT  TRANSFER  COEFFICIENTS  IN  OSCILLATING 
FLOW  THAT  SHOWED  ENHANCED  HEAT  TRANSFER 
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FIGURE  48  -  LOCAL  HEAT  TRANSFER  COEFFICIENTS  IN  OSCILLATING 
FLOW  THAT  SHOWED  DEGRADED  HEAT  TRANSFER 
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FIGURE  49  -  LOCAL  HEAT  TRANSFER  COEFFICIENTS  IN  OSCILLATING 
FLOW  THAT  SHOWED  DEGRADED  HEAT  TRANSFER 
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FIGURE  50  -  LOCAL  HEAT  TRANSFER  COEFFICIENTS  IN  OSCILLATING 
FLOW  THAT  SHOWED  DEGRADED  HEAT  TRANSFER 
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FIGURE  51  -  LOCAL  HEAT  TRANSFER  COEFFICIENTS  IN  OSCILLATING 
FLOW  THAT  SHOWED  DEGRADED  HEAT  TRANSFER 
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TABLE  I 
SUMMARY  OF  AVERAGE  HEAT  TRANSFER  RESULTS 


Run 

Yaw 

Average 

umber 

Angle 

Nu 

NRe 

NF 

Ns 

NA 

5 

90° 

565.49 

314057 

0.0 

0.0 

0.0 

6 

90° 

380.50 

157028 

0.0 

0.0 

0.0 

7 

90° 

304.32 

109598 

1.736 

X 

io"6 

0.190 

0.105 

8 

90° 

431.65 

213959 

0.464 

X 

lO"6 

0.099 

0.180 

9 

90° 

517.54 

283129 

1.166 

X 

IO"6 

0.330 

0.115 

10 

90° 

359.58 

151994 

3.962 

X 

io"6 

0.  602 

0.060 

11 

90° 

415. 10 

154470 

8.859 

X 

IO"6 

1.369 

0.150 

12 

90° 

387.19 

154470 

17 .719 

X 

IO"6 

2.737 

0.070 

13 

90° 

309.30 

149596 

0.914 

X 

io"6 

0.137 

0.100 

15 

90° 

530.16 

204846 

9.861 

X 

IO"6 

2.020 

0.052 

16 

90° 

418.60 

149596 

23.053 

X 

IO"6 

3.449 

0.340 

17 

90° 

594.91 

232762 

9.725 

X 

IO"6 

2.264 

0.690 

18 

90° 

450.24 

212769 

1.115 

X 

IO"6 

0.237 

0.150 

19 

90° 

400.49 

152811 

2.149 

X 

IO"6 

0.328 

0.110 

20 

90° 

717.90 

269246 

2.885 

X 

IO"6 

0.777 

0.310 

21 

45° 

360.71 

155739 

3.866 

X 

IO"6 

0.852 

0.040 

22 

45° 

515.02 

254242 

1.439 

X 

IO"6 

0.517 

0.055 

23 

45° 

362.85 

154891 

17  .671 

X 

IO"6 

3.871 

0.045 

24 

45° 

460.74 

217174 

9.013 

X 

io-6 

2.768 

0.049 

25 

45° 

356.32 

156166 

0.876 

X 

IO-6 

0.194 

0.100 

26 

45° 

462.35 

221461 

0.436 

X 

IO"6 

0.137 

0.120 

27 

45° 

567.89 

314057 

0.0 

0.0 

0.0 

28 

45° 

343.13 

154052 

8.883 

X 

10~6 

1.935 

0.110 

29 

45° 

556.50 

280040 

2.685 

X 

io"6 

1.063 

0.  220 

30 

45° 

497.50 

259876 

8.030 

X 

IO'6 

2.951 

0.095 

31 

45° 

335.60 

155314 

21.852 

X 

IO"6 

4.800 

0.048 

32 

45° 

453.45 

203253 

1.023 

X 

IO"6 

0.295 

0.115 

33 

45° 

374.23 

155739 

1.757 

X 

IO"6 

0.381 

0.090 

34 

45° 

357.48 

155314 

0.0 

0.0 

0.0 
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TABLE    I    Continued 


Run 

Yaw 

Average 

Number 

Angle 

Nu 

NRe 

NF 

35 

45° 

289.59 

110894 

3 

.452 

X 

10 

36 

45° 

341.41 

149788 

3, 

.769 

X 

10 

37 

45° 

272.71 

107462 

1 

.846 

X 

10 

38 

45° 

256.67 

108953 

20, 

.024 

X 

10 

39 

45° 

273.36 

110647 

0, 

.0 

40 

45° 

122.07 

35630 

0, 

.0 

41 

90° 

419.91 

150387 

4, 

.004 

X 

10 

42 

90° 

350.94 

94381 

25, 

,149 

X 

10 

43 

90° 

348 .06 

105535 

3, 

,638 

X 

10 

N, 


N 


-6 
-6 
-6 


-6 

-6 

-6 


0, 

,541 

0, 

.440 

0, 

,  798 

0, 

.240 

0, 

,280 

0, 

.650 

3, 

,085 

0, 

.410 

0, 

,0 

0, 

.0 

0, 

,0 

0, 

.0 

0, 

,602 

0, 

.110 

2. 

,374 

0, 

.330 

0, 

,384 

0, 

.290 
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